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We solve the nonlinear Poisson-Boltzmann equation for two parallel and likely charged plates 
both inside a symmetric elecrolyte, and inside a 2 : 1 asymmetric electrolyte, in terms of Weierstrass 
elliptic functions. From these solutions we derive the functional relation between the surface charge 
density, the plate separation, and the pressure between plates. For the one plate problem, we obtain 
exact expressions for the electrostatic potential and for the renormalized surface charge density, 
both in symmetric and in asymmetric electrolytes. For the two plate problems, we obtain new exact 
asymptotic results in various regimes. 
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I. INTRODUCTION 

When a charged object is inserted into an electrolyte, 
it attracts ions of opposite charge and repels ions of like 
charge. This leads to the well-known phenomenon of 
screening: the total electrostatic potential, due to both 
the external charges and the electrolyte, is exponential 
damped as a function of distance from the charged object. 
This phenomenon was first studied by Debye and Hiickcl 
(DH) and is often called Debye screening. At the 
level of mean field theory, the screened potential satisfies 
the so-called Poisson-Boltzmann (PB) equation, which, 
for a symmetric electrolyte with ion charges ±q, and ion 
number density n for each specie, is given by 



eA(j) + 2nq sinh/3q<^ = 0, 



(1.1) 



where (3 is the Boltzmann factor, q = 1.6 x 10 _19 C the 
charge of an electron. It is convenient to use the dimen- 
sionless potential = f3q<j), and to measure length in 
terms of the Debye length 
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k- 1 = y/e/2Pqh 



q : -q. 



(1.2) 



If the surface charge density on the object is low so that 
| M/ 1 <C 1 everywhere, the PB Eq. (|1.1[) can be further 
linearized to yield: 



- A* + * = 0, 



(1.3) 



which admits the famous screened Coulomb potential 
e~ r /4Trr as its Green's function. 

Linearization is however not appropriate for strongly 
charged objects. On the other hand, the full PB equa- 
tion is difficult to solve due to its nonlinear nature. For 
the simple case of one positively charged plate of infinite 
size, the exact solution has been known since the time of 
Verwey and Overbeek (2[: 
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= 2Iog- 



(1.4) 



where z is the coordinate perpendicular to the plate. The 
potential by a negatively charged plate is just the nega- 
tive of Eq. (|1.4|) . Note that Eq. (|1.4[) has a logarithmic 
singularity at 2 = 0. 

The general solution to the one plate problem is 1 &(z + 
zq) with an arbitrary parameter Zq. Normally one would 
have to determine the constant zq with the surface charge 
density a and the position of the plate fixed. For our 
purpose, however, it is much more convenient to fix the 
electrostatic potential to be Eq. (jl.4p and let zq be the 
position of the plate, which should be adjusted according 
to the Neumann boundary condition 



tf'(zo) = PqtDB 



d(j){zo) 
dn 



(3qcr 



(1.5) 

Defining a Gouy-Chapman length Iqc an d a dimension- 
less surface charge density r\ via 
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e/q(3a 
f3qa = 

€K 



Gouy-Chapman length(1.6) 
(1.7) 



the boundary condition Eq. (|1.5[) can now be expressed 
into the following concise form: 



sinh zq 



(1.8) 



The solution consisting of Eqs. (|1.4|1.8|) . though simple 
enough, shows an interesting property. For an arbitrary 
given surface charge density i] > 0, Eq. (|1.8j) can always 
be solved for zq > 0. That is, there is a one-parameter 
family of systems that gives the same potential Eq. Ji.^[ ). 
Furthermore, we can let the surface charge density n ap- 
proach infinity, then Eq. (|1.8[) dictates zq — > 0. This 
means that Eq. is the electrostatic potential pro- 

duced by a plate with infinite surface charge density at 
the origin zq — 0. In Fig. [T] we illustrate two plates 
(with surface charge density determined by Eq. (|1.8|) ) 
that produce the same potential Eq. (|1.4p . 

In the near-field regime z C 1, the potential Eq. (|1.4[) 
reduces to the famous Gouy-Chapman solution: 
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*(z) = -21ogz + 0(l), near field. 



(1.9) 
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FIG. 1: Up: The potential Eq. (|1.4[1 can be produced by 
two plates with different surface charge densities and differ- 
ent locations. Both satisfy the Neumann boundary condition 
Eq. (|1.5p . The extremal possibility is an infinitely charged 
plate at zo — 0. Down: The renormalized surface charge den- 
sity r\R as a function of the bare surface charge density r\ on 
a single uniformly charged plate. Red solid line: a q : —q 
symmetric electrolyte as given by Eq. (|1.12p . r\R saturates at 
4 in both directions. Blue dashed line: a 2q : —q asymmet- 
ric electrolyte as given by Eq. (| 1 . ltjf) . r)a saturates at 6 and 
6(2 — \/3) in the positive and negative sides respectively. 



In the far-field regime z> 1, Eq. (|1.4[) can be expanded 
into Taylor series of e~ z : 



= Ae~ z +Q{e 



-2z\ 



4e" 



-Sz 



+ 0{e~ 2z ), (1.10) 



where Sz — z — z is the distance to the charged plate at 
zq. It is understood that the parameter zq is a function 
of 77 as determined by Eq. Q1.8p . On the other hand, in 
the linear PB theory, a plate with surface charge density 
t]r at the same location zq produces at z a potential 



linear 



(z) = r/ R e 



-Sz 



(1.11) 



Following the proposal by Alexander et. al. [|| , we iden- 
tify the far-field asymptotics Eq. (|1.10[) of the nonlin- 
ear theory with the linear theory Eq. (|1.11[) , and define 
the renormalized (or effective) charge density i]n. Using 



Eq. (|1.8[) to eliminate zq in favor of 77, we find the renor- 
malized surface charge density 77^ (77) as a function of the 
bare surface charge density 77: 



m(v) = 



2ri 



1 + V1 + W 



(1.12) 



which is illustrated by the red solid curve in the right 
panel of Fig. [TJ In the weakly charged limit 77 <C 1, the 
renormalized charge density 77^ is just the bare one 77^ — > 
77; in the strongly charged limit 77 3> 1, 77^ saturates at 
t]h(+(X)) = 4. r/ R saturates at —4 when 77 is negative and 
large. More generally, Eq. (|1.12[) is invariant under the 
charge-inversion transformation (77, r)n) — > (— 77, — T] R ), in 
accordance with the fact that the electrolyte is symmet- 
ric. The practice of using linear theory with a renor- 
malized, i.e. effective, surface charge density in the far 
field is usually called charge renormalization following the 
seminal work by Alexander et al. Here the most strik- 
ing property is that 77^ saturates at a finite value in the 
strongly charged limit. It shows that electrolytes are able 
to screen infinitely charged objects within finite distance. 
This should be regarded as one fundamental property of 
the Poisson-Boltzmann theory. 

In general, the renormalized surface charge density de- 
pends on the shape of the charged object. It also depends 
on the properties of electrolyte. For example, our anal- 
yses in this work show that inside a 2q : —q asymmetric 
electrolyte, a positively charged plate produces a dimen- 
sionless potential 



2q:-q 



(z) = log 



l + 4e- 



-2.1 



{l-e- z f 



(1.13) 



which diverges logarithmically at z = 0. A negatively 
charged plate in a 2q : —q electrolyte, on the other hand, 
produces a potential 



2<i 



:- g 0) = log 



l-4e" 



(1 



-z\2 



(1.14) 



which diverges at z m = log (2 + \/3). These potentials 
were discovered by Andrietti et al in 1976, but has re- 
mained largely unknown since then. The far-field asymp- 
totics of these potentials are given by 



2<l 



: _ q (z)~±6e- z + 0(e~ 2z ). 



(1.15) 



Carrying out a similar analysis as in the case of sym- 
metric electrolyte, we find that the relation between the 
renormalized surface charge density 77^ and the bare sur- 
face charge density 77 as 



36r] R (n R + 6) 



(6 



m) (vr 



2A m + 36) 



(1.16) 



which is illustrated by the dashed blue curve in Fig. [T] 
It is interesting to note that for 77 > 0, ?tr saturates at 
t]r(oo) = 6, while for 77 < 0, r\ R saturates at 77^— 00) = 
-6(2 - s/3) w -1.6077. The fact 77^(77) ^ -77^-77) 
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FIG. 2: The two dimensional parameter space for the two 
plate problem. There are four regimes where simple asymp- 
totic results can be obtained without solving the nonlinear 
PB. Our solution in terms of Weierstrass elliptic functions re- 
duces to the known asymptotic results in the corresponding 
regimes, and works in the cross-over regions shown shaded in 
the figure as well. Part of this plot was taken from the review 
by Andelman 0). The blue line: two strongly charged plates 
with a variable distance. The red line: two weakly charged 
plates with a variable distance. 

shows that the electrolyte is indeed asymmetric. There 
does not seem to exist simple analytic result for the one 
plate problem in more general mq : —nq asymmetric elec- 
trolytes. 

The finiteness of electrostatic potential away from an 
infinitely charged object in electrolyte is a short scale 
property of the PB equations, and therefore is expected 
to hold independent of the large scale geometry of the 
charged objects. Indeed the phenomenon of charge renor- 
malization was first discovered in numerical studies of 
spherical colloids It was found there that the effec- 
tive interaction between charged objects appears to re- 
main finite when the surface charge density is scaled up. 
However that has not been any analytic result showing 
that this is indeed true [l3| . For the case of spherical col- 
loids, even the potential by a single colloid is unknown. 
In general, lack of analytic results is one major obstacle 
to understanding of the physics of charged colloids. 

The problem of two charged plates is the simplest toy 
model for effective interactions between charged colloids. 
For two likely charged inside a symmetric electrolyte, a 
formal solution of the nonlinear PB was obtained in terms 
of elliptic integrals in the classical monograph by Verwey 
and Overbeek Q as early as 1940's. From this, asymp- 
totic results can be obtained in four regimes shown in 
Fig. [2j For a detailed discussion, see the review by An- 
delman [i| . There however has been no analytic result in 
the cross-over regions, shaded in Fig. [2J In the present 
work, we shall analytically calculate the electrostatic po- 



tential profile and the interaction between the two iden- 
tically charged plates, both in symmetric and in 2 : I 
asymmetric electrolytes. By expressing the electrostatic 
potential in terms of Weierstrass elliptic functions, we 
are able to calculate the interaction between two plates 
in the whole parameter space. We derive exact results in 
various asymptotic regimes as well as cross-over regimes. 
In particular, we find analytic expressions for the interac- 
tions between two infinitely charged plates, both in sym- 
metric and asymmetric electrolytes. These results consti- 
tute the most direct demonstration of charge renormal- 
ization. We also find that the potential by a single plate 
in a 2q : —q asymmetric electrolyte in terms of elemen- 
tary functions, as discussed above. Last but not least, 
our result also demonstrate how bi-valence counter-ions 
dramatically reduces the interaction between two likely 
charged plates. It is important to note that our solution 
should be equivalent to one by Verwey and Overbeek [2[ . 
In fact, it is known that all elliptic functions can be ex- 
pressed in terms of Weierstrass elliptic functions. The 
representation in terms of Weierstrass elliptic functions 
is however particularly convenient because of their well 
known analytic structure. It allows us to deduce many 
analytical properties of the PB equation for two-plates 
geometry. 

In recent years, there has been much discussion on the 
physics of charged colloids beyond Poisson-Boltzmann 
theory Of particular interests are the counter- 

intuitive possibilities such as charge inversion, like-charge 
attraction, etc. Partly due to the lack of reliable analytic 
results, however, the situation remains murky after more 
than a decade of debates. In this work, we shall limit 
ourself inside the domain of Poisson-Boltzmann theory. 
It is our belief that the physics beyond PB can be un- 
derstood only after we understand the physics within PB 
theory. 

The remaining of this paper is organized as follows. In 
Sec. [n] we construct the solution to PB in terms of Weier- 
strass elliptic function, both in symmetric electrolytes 
and in 2q : — q asymmetric electrolytes. We also ex- 
press the interaction between two plates in terms of the 
first integral of the PB equation. In Sec. IIII1 we discuss 
the one plate problems, analytically calculate the elec- 
trostatic potential and the renormalized surface charge 
density as a function of the bare surface charge density. 
In Sec. El we study the problem of two infinitely charged 
plates. We obtain the exact asymptotics both in the large 
separation regime and in the small separation regime, 
and find that multivalent counter-ions significantly re- 
duce the repulsion between two likely charged plates. We 
also obtain analytic results for arbitrary separation. In 
Sec. IVII we study two weakly charged plates and obtain 
exact asymptotics both in the large separation limit and 
in the small separation limit. In Sec. I VIII we analyze the 
case of two strongly charged plates and briefly discuss 
the most general case where the surface charge density is 
neither strong nor weak. Finally in Sec. I Villi we summa- 
rize our result and discuss possible future directions. In 
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Sec. [A] we present two different derivations of Poisson- 
Boltzmann equation in symmetric electrolytes, and dis- 
cuss some subtleties in each derivation. In Sec. [5] we 
discuss some general properties of PB equation with the 
geometry of two parallel charged plates. The pressure 
between plates naturally emerges as the first integral of 
PB equation. In Sec. [C] we derive PB for the case of 
asymmetric electrolytes and discuss some of its general 
properties. 



II. SOLVING POISSON-BOLTZMANN 
EQUATION 

Inside a symmetric electrolyte, the Poisson-Boltzmann 
equation for the dimcnsionless electrostatic potential = 
between two identical charged plates is given by 



+ sinh*(z) = 0. 



(2.1) 



Let e, e' be the dielectric constants for the electrolyte and 
for the plates. In the limit e>e', the relevant electro- 
static boundary condition reduces to that of Neumann: 



9* 

dn 



plates 



where the normal director points into the plates, 
easy to obtain the following first integral: 



(2.2) 



It is 



a = cosh*(z) - i(*'(z)) 2 



cosh tf(MP), (2.3) 



where "MP" stands for "middle point" . The net interac- 
tion between two plates only depends on the constant a 
and is given by 



P not = 2nT(a-l) 



TSa 



4ire 2 DB \ Bj 



(2.4) 



where 8 a = a — 1, and 2n is the total density of ions in 
the bulk, a > 1 for two likely charged plates, a — >• 1 as 
the distance between plates becomes large. For detailed 
derivations of all these results, see Sec. lAl 

Now introducing a new function p{z) through 



4p(z) + 2a/3 = e* (z) , 
the first integral Eq. (12.3[) is transformed into 
(p') 2 =4p 3 -g 2 p-g 3 , 

where 



<I2 



a 



i 

4' 



.93 = 



27 



Q 

24' 



(2.5) 



(2.6) 



(2.7) 



Eq. (|2.6[) is the differential equation satisfied by the 
doubly periodic Weierstrass elliptic function p{z\ g 2 , 53), 
with g 2 and g 3 two of its invariants. Below we invoke 
some of its well known properties without any proof. The 



readers are referred to the classic monograph by Whit- 
taker and Waston Q for relevant details. The function 
p{z;g 2 ,g 3 ) can be explicitly represented as a double se- 
ries in z: 



1 



1 



(z — 2mu)i — 2nuj 2 ) 2 
1 



(2mu)i + 2nuj 2 y 



(2.8) 



where the prime in the summation means exclusion of the 
term with m = n = 0. p(z) is a mcromorphic function of 
z, generally treated as a complex variable. 2uj\ and 2^2 
(complex numbers in general) are the two periods and 
are related to the two invariants via 

. V 1 

g 2 (oj!,uj 2 ) = 60} t 

m.n 



(2mcji + 2nuj 2 ) 4 



g 3 (uj u uj 2 ) = 140 



1 



(2mwi + 2nco 2 ) 6 



(2.9a) 
(2.9b) 



All series involved here are absolutely convergent. The 
function p(z; g 2l g 3 ) is completely determined by two in- 
variants 32,33, or cquivalcntly, its two periods 2lui,2lu 2 . 
In our case, they depend on only one parameter a, see 
Eq. (|2.7p . One can easily check using Eq. (|2.8[) that p{z) 
is indeed periodic with respect to two periods 2u>\ and 
2w 2 : 



p(z) = p(z + 2wi) = p(z + 2u} 2 ). 



(2.10) 



Furthermore, two periods are guaranteed to be linearly 
independent as two dimensional vectors in the complex 
plane @. From Eq. (|2.8[) . z = is a second order pole 
of p(z). By periodicity, p[z) has an infinite number of 
second order poles z m ^ n = 2muji + 2nLtj 2 , which form a 2D 
oblique lattice. These are the only singularities of p{z). 
In the remaining of this work, we shall focus on the case 
of two likely charged plates where a > 1. For this case, 
it can be explicitly shown that the lattice of singularities 
is rectangle. Hence 2wi and 2uj 2 can be chosen to be 
real and purely imaginary respectively. This lattice is 
illustrated in Fig. [3] The case with a < 1 (two oppositely 
charged plates) is slightly more complicated and will not 
be discussed here. For a > 1, the real half period u>i is 
related to a through the two invariants Eq. (|2.7[) : 



wi(a) 



Twhp (52 (a), 53(a)) 
« 2 1 a 3 



WHP 



£L_i,£L_fLj, (2.11) 

3 4' 27 24.'' v ; 



where the subscript "WHP" stands for "Weierstrass Half 
period" . This function can be conveniently computed 
using Wolfram Mathematica 7. [l4| 

From now on we shall consider real variable, 

since this is ultimately what we concern about. Inverting 
Eq. (|2.5j) . we find that the dimensionless potential ^ is 
given by 



#(z) = log 



1 a" 
4' 27 



2a 



(2.12) 
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FIG. 3: (Up) The periodic lattice structure of the singularities 
of p(z, gi, gs) in the complex plane. Each grid point is a 
pole of second order. For a > 1, 2cji is real, while 2cj2 is 
purely imaginary. (Down) Black solid curve: Eq. (|2.12p as 
a function of normalized real coordinate z/2wi with a > 1. 
2 = 0, 2ui\ are logarithmic singularities; z = Wi is mid point 
between plates. Here a = 10, corresponding to a period 2oji ~ 
1.40761 Idb- Two dashed lines are two maximally separated 
plates with infinite surface charge density, for given pressure 
a. Two solid (red) lines are two plates with finite surface 
charge density that gives the same pressure. The parameter 
zo in Eq. (|7.1I) is the distance between the dashed line and 
the solid line, i.e. it is the distance between the plate and 
the logarithmic singularity nearby. Blue dashed curve: the 
potential Eq. (|2.21[) inside 2q : —q asymmetric electrolytes 
with q = 1, and 2wi — 4.52566 Idb- 



densities in such a way that the Neumann BC is always 
satisfied and the constant a is fixed. That is, there is 
a one parameter family of problems that the yield the 
identical potential profile. Since the potential between 
two charged plates must be a smooth function, the two 
singularities 0,2wi, must be outside two plates. This 
means that for a given pressure a, there is a maximum 
separation between two plates, which is precisely the 
real period 2ui\ . The surface charge density yielding this 
interaction at this maximal separation is infinity, see the 
dashed vertical lines in Fig. [3J 

Similar to the one plate problem, we let the plate on 
the left to be at zq. Since the other plate carries the 
same charge density, and since p(z) is an even function, 
the other plate must be located at z = 2u>\ — zq . The val- 
ues of a, zo should then be determined by two boundary 
conditions Eq. (|2.2[) . 

By our choice of coordinate system, z = u>\ is the mid 
point between two plates. It is known that p at this point 
satisfies the following equation Q: 



At 6 - g 2 t - g 3 = 0, 



(2-13) 



with <72 j 53 given by Eq. (|2.7p . This equation admits three 
real roots if a > 1: 



1 



(2.14) 

Only the last one gives a positive potential "J through 
Eq. (|2.12j) . It is also easy to check that this solution 
satisfies Eq. ((231) . 

We note in passing that from the electrostatic potential 
Eq. (|2.12j) . one can straightforwardly calculate the den- 
sity of counter-ions and co-ions between the plates. We 
shall however not elaborate on this issue in the present 
work. 



A. 2q : —q Asymmetric Electrolyte 



which has logarithmic singularities at z = 
0, ±2wi, ±4wi, . . .. These singularities have the same 
nature as the one in the solution for the one-plate 
problem. In fact, Eq. (| 1 .4[) is an elliptic function with 
(2aJi,2uj 2 ) = (+oo,27ri). In Fig. [31 we plot this potential 
in one period along the real axis z <G (0, 2wi) for a 
particular value of a = 10. It is positive everywhere, 
and diverges logarithmically at z = 0, 2lj\ w 1.4076. It 
therefore describes the potential between two positively 
charged plates. The potential between two negatively 
charged plates is simply the negative of Eq. (|2.12[) . since 
the saddle point equation Eq. (]2.ip is invariant under 
transformation , J7 — > — "t. This symmetry no longer 
holds in an asymmetric electrolyte, as we will see shortly 
after. One possible loci for two positively charged plates 
arc illustrated by the solid lines in Eq. ([3]). Like the 
case of a single charged plate, we can continuously 
change the loci of the plates as well as the surface charge 



Let positive ions and negative ions carry charge of 2q 
and — q respectively To keep overcall charge neutrality, 
their bulk ion densities are n/2 and n respectively. The 
Debyc length for a 2q : ~q electrolyte is [l(| 




(2.15) 



The Poisson-Boltzmann equation in a 2q : —q electrolyte 
is 

-d^+ 1 -e*- 1 -e-™=Q. (2.16) 

Linearization of this equation leads to Eq. (|1.3p . as it 
should be. The first integral of Eq. ([2~T6]) is: 

= L •(«*> + l e -mMP) (217) 
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For likely charged plates, a > 1/2. a converges to 1/2 
as the distance between plates becomes large. The net 
interaction between two plates is given by 



Pu 



a-rl 



TSa 



(2.18) 



where 5a = a — 1/2. For detailed derivation of the above 
results, see Sec. [C] 

Introducing a variable p(z) via 



6 p(z) + a = cxp <I>, 



(2.19) 



we find that Eq. (|2.17[) reduces to Eq. (|2.6[) with two 
invariants given by 



92 



a 

Y 



93 



27 



1 

108' 



(2.20) 



Inverting Eq. (|2.19p , we find that the electrostatic poten- 
tial is given by 



^ = log 



6p (z 



a 2 a 3 
T' 27 



1 

108 



(2.21) 



For real z and a > 1/2, this potential is always positive. 
In Fig. |3]we plot Eq. (|2.21[) within one period on the real 
axis for the particular case of a = 1. Therefore Eq. k2.21\) 
with z on the real axis describes the electrostatic potential 
between two positively charged plates, as long as a > 1/2. 
The real period 2wi as a function of a is given by 



2wi(a) = 2 T WH p 



a 2 a 3 
~3~' 27 



1 

108 



(2.22) 



Because the electrolyte is asymmetric, the saddle point 
equation Eq. (|2.16[) is not invariant under the transfor- 
mation ^ — > — The potential by negatively charged 
plates has to be studied separately. We shall come back 
to this problem below. 



III. ONE PLATE PROBLEM 

For the case of symmetric electrolyte, a — > 1 as the 
distance between plates approaches infinity, see Eqs. (|2.3[ 
The two invariants Eq. (|2.7[) reduce to 



1 

92 -> 92 



1 

216' 



(3.1) 



Using Wolfram Mathematica 7, we can explicitly check 
1w\ — > +oo, 2ll> 2 — > 27TI. (3.2 
The electrostatic potential Eq. (|2. 12p then becomes 

1 1 \ 2 



q:-q 



(z) 



log 



12' 216 



(3.3) 



At a given finite z, this potential is completely due to 
the plate near z = 0, since the other plate is infinitely 



far away. On the other hand, we already know that the 
potential by one charged plate inside a symmetric elec- 
trolyte is given by Eq. (| 1 .4[) . Comparing it with Eq. 
we find the following useful identity: 



P z 



1 



1 



1 (1 



12' 216/ 4 \l-e 



(3.4) 



A. + plate in 2q : —q electrolyte 

Now consider a 2q : — q asymmetric electrolyte. When 
two plates are widely separated, a — > 1/2, sec Eqs. (|2.171 
12.181) . Interestingly, the two invariants Eqs. (|2.20[) ap- 
proach the same limit Eq. (|3.1|) . Using Eq. (|3.4[) and 
Eq. (|2.21l) . we find that inside a 2q : — q asymmetric elec- 
trolyte, the electrostatic potential of a positively charged 
plate is given by 



*2 q :- q (z) = log 



= log 



6p 



1 + 4er 



12' 



216 J 

,-2x 



(l-e~ z ) 2 



(3.5) 



which we have already presented in Eq. (|1.13[) . It is 
easy to check explicitly that Eq. (|3.5p indeed satisfies PB 
Eq. (|2.16l) . Furthermore for z > 0, Eq. (|3.5|) is positive 
everywhere. 

Let the plate be at zq with a surface charge density 
rj > 0, the boundary condition Eq. (|2.2[) dictates 



dz H H 



3(e 20 + 1) 



(e z ° - 1) (coshzo + 2) 



V . (3.6) 



On the other hand, the far field asymptotics (z 1) of 
Eq. (£03) is 



^2q:- q ( z ) ~ 6e_Z = 6e" z °e" fc , Sz = z-z . (3.7) 



Comparing this with the lienarized PB theory Eq. (|l.lll) , 
we find that the renormalized surface charge density rjR 
is given by 



7} R = 6er Zo > 0, 



(3.8) 



where z is in turn related to rj by Eq. (|3.6p . In the limit 
?; — > oo, zo — > + according to Eq. p.6[) , and 77^ saturates 
at 6. Therefore Eq. p.5[) ^iwes </ie potential by a positive 
infinitely charged plate located at z — 0. For finite value 
of 77, we can use Eq. (|3.8p to eliminate zo in favor of 77 in 
Eq. (|3.6p . and obtain the relation between 77^ and 77: 



36 r] R (r] R + 6) 



(6 -VR) (^ + 2477/? + 36) 



= V, 



(3.9) 



which we already presented in Eq. (|1.16p in Sec. |TJ It is 
easy to verify using this relation that 77^ — > 6 as 77 — > 00. 
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B. — plate in 2q : — q electrolyte 

Rather than presenting the derivation that leads to the 
correct solution, we directly show the result here. The 
function 



#o„._„(z) = ^t.^Jz + in) 



2q:-q 



log 



q:-q\ 

1 - ie~ z + e 



-2z 



(1 



(3.10) 



also satisfies the differential equation Eq. (|2.16[) . and is 
negative as long as z > log (2 + \/3) rj 1.317. It has 
a logarithmic singularity at z = log(2 + v3). It there- 
fore describes the potential by a negatively charged plate 
inside a 2q : — q electrolyte. 

Let the plate be at zq with a surface charge density 
i] < 0, the boundary condition Eq. (|2.2[) dictates 



(*o) 



3(e 2 ° - 1) 



(e z ° + 1) (cosh;z 



r,. (3.11) 



In the limit r\ — > — oo, we have Zo — > cosh -1 2 = log (2 + 
V3). Therefore Eq. p.lOp describes the potential by a 
plate with r\ = — oo at zq = log(2 + \/3). In the far 
field z 3> 1, the asymptotic behavior of the potential 
Eq. (jBTDl is 



(z) rj 6e~ 2 = -6e" 



z - (3.12) 



Comparing this with the lienarized PB theory Eq. (jl.lll) , 
we find that the renormalized surface charge density rjn 
is given by 



-6e" 



< 0, 



(3.13) 



where zq is in turn related to r\ by Eq. (|3.1ip . Eliminating 
Zo in favor of r\ in Eq. (|3.11[) using Eq. (|3.13[) . we find the 
relation between the renormalized surface charge density 
r]n and the bare density i] (both negative in this case), 
which is identical to Eq. (|3.9[) . rjn as a function of r\ is 
illustrated in Fig. [T] for both positive and negative rj. 
T]r saturates at —6/(2 + \f?>) sa 1.6077 as r] — >• — oo, and 
saturates at 6 as rj +oo. 



IV. TWO WIDELY SEPARATED PLATES WITH 
ARBITRARY CHARGES 

The solution to the one plate problem help us under- 
standing the asymptotics of the two plate problem where 
the separation is much longer than the Debye length, 
L 3> £db- The potential in the middle is approximately 
the sum of that of two isolated plates. In the middle of 
plates, therefore, 



(4.1) 



Using Eqs. (gHH Ei, or E( l s - <EH3 I2TT5]) . and the fact 
that ^(MP) <C 1, the net interaction between two plates 



is approximately 
da 

"net 



1 -^{MP) 2 = 2 m (rj) 2 e- L /^ 



2ir£ 2 DB \ Bj 



(4.2) 



where we have restored the physical unit for length, while 
r/n is given by Eq. (|1.12[) for symmetric electrolytes and 
by Eq. (|1.16[) for 2q : —q asymmetric electrolytes. Note 
that this result be obtained from the corresponding result 
Eq. (|6.11l) in the Debye-Hiickcl regime by replacing the 
bare surface charge density rj with the renormalized one 
r)R in the regime where L £jjb ■ This result applies to 
the whole region to the left of the vertical line £db /L = 1 
in Fig. [21 which includes the intermediate regime as well 
as a major part of the Debye-Hiickel regime. 



V. TWO INFINITELY CHARGED PLATES 

A. Large Separation Asymptotics: "Intermediate 
Regime" 



Taking the limit r] — > ±oo in Eq. (|4.2[) and using 
Eqs. (|1.12[ 13. 9p . we obtain the large separation asymp- 
totics for the interaction between two infinitely charged 
plates: 

ST „-L/£ DB 



18T 
— vi — ^ 1 



-L/t DE 



±00 plates in q : —q. 
+oo plates 2q : —q; (5.1) 



18(2-y/3) 2 T Lllnr 



— ^ — e ~i-~ljb plates 2q : —q. 

Note the asymmetry between positively and negatively 
charged plates in asymmetric electrolytes. 



B. Small Separation Asymptotics: 
"Gouy-Chapman Regime" 

When L <C Idb, the potential at the middle is high, 
so that the co-ion density is negligibly small between two 
plates. The corresponding term is exponentially small 
and can be dropped from the PB equation. The first 
integral then approximately reduces to 

ie*-i(9 2 *) 2 , ±oo plates q : -q 



\{d z ^) 2 , +oo plates 2q : -q , (5.2) 



6 C 2 1 



\{d z ^) 2 , — oo plates 2q : —q 
which can be integrated easily to yield: 

2 lo s tMci> 00 P lates i '■ -<i 

*(z) = { 2 log -^fr + log §, +oo plates 2q : -q , 



-log 



3S Kz 

4K 
os 2Kz 



\ log |, — oo plates 2q : —q 



(5.3) 
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L 




FIG. 4: Log-linear plots of L v.s. 8a for two infinitely charged 
plates. Black: q : —q symmetric electrolyte, as given by 
Eq. (|5.6|) . Red: r\ = +oo, 2q : —q electrolyte, see Eq. (|5.7|) . 
Green: rj = -co, 2q : —q electrolyte, see Eq. (15.91) . 



6 P(z+ai 2 )+a 
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where A' = y/u/2. vO'(z) diverges at z ~ ±n/2K for the 
first two cases, and at zLtt/AK for the last case, therefore 
the distance betwee n two infinitely charged plates is given 
by L = n/K = n^2/a for the first two cases and L = 
7r/\/2a for the last case, a as a function of L is then 
given by 



2ir 2 /L 2 , ±oo plates q : —q 
2tt 2 /L 2 , +oo plates 2q : —q 
7r 2 /2L 2 , — oo plates 2q : —q 



(5.4) 



Using Eq. (|2.4[) and Eq. (|2 . 18|) . the pressure, in physical 
unit, is then given by 



P 



2-rr 2 eT 

2ir 2 eT 

7r 2 eT 
2<r ! L 2 ' 



±co plates q : — q 
+oo plates 2q : — q 
— oo plates 2q : —q 



(5.5) 



which is independent of the co-ions and the surface charge 
density. 



C. Arbitrary Separation 



Inside a symmetric electrolyte, Eq. (|2.12[) gives the po- 
tential between two infinitely charged plates located at 
z = 0, 2<ji, with uj\ given by Eq. (|2.11j) . Replacing ui\ 
with L/2 in Eq. (|2.1ip . we find the relation between L 
and a: 

fa 2 la 3 a \ 
L = 2 Twhp I y - 4> 27 - 24 ) ' ±0 ° p es q : 

(5.6) 



FIG. 5: Up: The function 6p(z + u)2\g2,gz) + a has two 
zeros (Ci,C2) within one period. Between these two zeros, 
< 6p + a < 1. For the case shown in this plot, a — 
5/9, (2wi,2w a ) = (7.084, 6.225i), (Ci,Ca) = (1.275,5.810). 
Down: For £i < jz < * = log(6 p + a) < describes 
the potential between two negatively charged plates. The 
potential has two logarithmic singularities at (£1,^2). 



Inside a 2q : — q asymmetric electrolyte, the potential 
Eq. (|2.21[) is always positive and has logarithmic singu- 
larities at z = 0, 2wi as long as a > 1/2. It therefore 
gives the potential between two r/ = +00 plates located 
at 0, 2wi(a), with two invariants given by Eq. (|2.20p . The 
relation between L and a is then given by 



L = 2T 



WHP 



a 2 a 3 
Y' 27 



1 

108 



-00 plates 2q : —q. 



(5.7) 

The potential between two negatively charged plates 
inside a 2q : —q asymmetric electrolyte is also given by 
Eq. (|2.21[) . but with z not lying on the real axis. In fact, 
we find that for a > 1/2, 6 p(z + W2; 52,33) + a (with 
z a real variable) has two zeros £1, (2 within one period. 
Between these two zeros, < 6p + a < 1. Note that 
UJ2 is always purely imaginary for a > 1/2. In Fig. [5] we 
illustrate the function 6 p(z+w 2 ; 32, ff3)+Q : with a = 5/9, 
(2wi,2w 2 ) = (7.084,6.225z).The function has two zeros 
(C11C2) = (1-275,5.810) within one period. Since p is 
an even function of z, we have £2 = 2u\ — Ci- As is 
also shown in Fig. [5l the potential $ = log(6 p + a) 
has two logarithmic singularities at (£1, £2); it is negative 
between these two singuarlties. vff therefore describes the 
potential between two r\ = —00 plates located at Ci>Ca 
respectively. For the case shown in Fig. [SJ the distance 
between two plates is (2 — Ci = 2wi — 2£i = 4.535. 

For arbitrary a, there is no analytic expression for the 
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two zeros (Ci, (2)- Nevertheless, we can alway find them 
numerically using Mathematica by solving the following 
equations: 

/ a 2 a 3 1 \ 

6pl{i(a)+u 2 (ay,—, — -— j +a = 0,(5.8) 

< Ci(a) < (2(a) < 2wi(a). 

The relation between L and a for a pair of negative in- 
finitely charged plates is given by 

(a 2 a 3 1 \ 

—00 plates 2q : —q. (5-9) 
The relations Eqs. (EH EH ED arc plotted in Fig. g] 



VI. TWO WEAKLY CHARGED PLATES 

A. Small Separation Asymptotics: "Ideal Gas 
Regime" and Beyond 

Here we focus on the case of symmetric electrolytes 
only. If L < 1, and 77 < 1, the potential Eq. ([232]) 
between plates can be expanded into Taylor series around 
z = cui up to Sz 2 : 

V(z) = *(wi) + ^"(w^Sz 2 + 0{5z A ). (6.1) 

The condition under which the expansion is quantita- 
tively accurate will be clear below, ^"(cji) can be calcu- 
lated using Eqs. (j2~T2l) . Eq. ([22]), and Eq. (f2~T4l) : 



= \Ja 2 - 1. 



(6.2) 



The Neumann boundary condition Eq. (12.21) (evaluated 
at 8z = L/2) then gives: 



a = \ 1 



4?7 2 
~p7 



<7 : -9- 



(6.3) 



The net interaction between plates can be calculated us- 
ing Eq. (J2U): 



4n£ 2 DB \ Bj 



\iacLj 



q : -q. 



(6.4) 

The Taylor expansion works if ty"(wi)L/2 -C 1. Using 
the solution Eq. (|6.3|) we find that this condition trans- 
lates into r\L <C 1, or in physical units, 



L < £gcJdb, 



(6.5) 



which is the top right quarter in Fig. [2] This result has 
not been derived previously. 

This quarter can be divided into two regimes: if 
l 2 DB / IgcL ^> 1, we have a w 2rj/L, which is the ideal 



gas regime [J], shown in Fig. [2j The net interaction in 
physical units is then approximately given by 



P 



T 



2Ta 



2nX£ GC L qL 



(6.6) 



In this regime, all co-ions are excluded from the region 
between plates by the strong Coulomb repulsion. The 
counter-ions, on the other hand, behaves as an ideal gas 
with number density 2a /(qL). If l 2 DB j IqcP "C 1, which 
is part of the Debye-Hiickel regime where the separation 
is much shorter than the Debye length L -C (db, 5a rj 
2rj 2 / L 2 . The net interaction in physical unit is then give 
by 



TP 2 
1 Z DB 



2tt\ 



Bj l GC 



L 2 ' 



(6.7) 



Both sub-regimes are discussed in the review by Andel- 
man 

Similar analyses can also be carried out for the cases 
of asymmetric electrolytes. The results however are 
rather complicated, and not particularly illuminating. 
We therefore only present asymptotic behaviors in the 
ideal gas regime, where L <C lac an d £ 2 db /£gcL ^S> 1: 



2To 
qL 



-plates 2q 



-q, 



^jj —plates 2q : —q. 



(6.8) 



The asymptotics of potential in the Debye-Hiickel regime 
will be discussed below. 



B. Large Separation Asymptotics: "Debye-Hiickel" 
Regime 

Now consider two weakly charged plates separated by 
a distance that is not too small, so that the potential is 
small everywhere f C 1. The Poisson-Boltzmann equa- 
tion can then be linearized. This is the Debye-Hiickel 
regime, where nonlinearities become irrelevant and dif- 
ferent electrolytes yield the same interaction and poten- 
tial profile, as long as they have the same Debye length 
and temperature. The linearized PB can be easily solved 
to yield the electrostatic potential: 



coshz 
sinhi/2' 



(6.9) 



where two plates are located at ±P/2 respectively, a is 
given by 



5a = 



2sinli (P/2) 



(6.10) 



Using Eq. (|2.4[) and Eq. ()2.18[) . the net interaction be- 
tween two plates is 



Tr, 2 

P not = , 2 7 smh 

0~Kt DB ABj 



2f 



DB 



.11) 
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The far field asymptotics is 
T 



Pne.t — 



2-K\ B jt GC 



DB- 



(6.12) 



Note that Eq. (g^J) can be obtained from Eq. (pTL2"j) by 
replacing r\ with r)R. 



VII. TWO STRONGLY CHARGED PLATES, 
AND THE GENERAL CASE 



In the general case of finite surface charge density and 
arbitrary separation between the plates, we need to de- 
termine two parameters a, zq, for a given 77 and L. For 
this purpose we need two equations. Note that zq by 
definition is the distance from the plate to the nearest 
singularity of p(z) z — 0, see the plot of in Fig. [3] 
Since two plates are located at zq, 2lu\ — zq and are sep- 
arated by L, we obtain the first equation: 



L = 2wi(a) - 2z , 



(7.1) 



where wi (a) is given by Eq. (|2.11[) for the case of symmet- 
ric electrolyte and by Eq. (|2.22j) for the case of 2q : —q 
asymmetric electrolyte. The other equation is obtained 
using the Neumann boundary condition Eq. (|2.2I) , which 
can be transformed into 



-p'(^o) = V p( z o) + 



■tx> 2 ) = V 



p{zo 



-q, 



. a 
; 6 



or ?/ > 2q : —q, 
(7.2a) 

, 7] < 2q : -q, 
(7.2b) 



after using of Eq. ff2TT2|) or Eq. ([231) . Note that both p 
and p' depend on a through two invariants 32(0;), 53(a), 
see Eq. (|2.7[) or Eq. (|2.20p . while the prime in p' is with 
respect to z. We are unable to obtain analytic solution to 
Eqs. (|7.1l 17.21) in the general case. Nevertheless, simple 
result can be obtained if 77 1. 

Let us first consider the first two cases given by 
Eq. (|7.2a[) . Since z = for 7/ = 00 (sec Fig. [3]), we 
expect zo <C 1 for 77 ^ 1 . Near zq = 0, the elliptic 
function in Eq. (|7.2p can be expanded in terms of zq as 



p(z ) = -s+0(«o)- 



(7.3) 



Note that the leading order term is independent of a. 
Hence Eq. (|7.2a[) reduces to 



2 1 1 

4 4 6 

r 1 =-(l + z 2 a/6)- 1 ^ 

ZQ 



2 

zo' 



Plugging it back into Eq. (|7.1[) we find that ?/ ^> 1, 

4 



2wi(a) = L + 2z ~ i 



g : —q, or 77 > 2q : — g. 

(7.4) 



Now consider Eq. (|7.2b[) . Since the right hand side van- 
ishes linearly at z = Ci , see the discussion after Eq. (|5.7|l , 
we can expand it around £i : 



(zq + uji) + -a = a (zo - Ci) + 



Substituting this back into Eq. (|7.2bl) we find 

1 



(7.5) 



(7.6) 



2wi(a) = L + 2Ci--, Ti<0 2q:-q. (7.7) 



z ~ Cl 



Combining this with Eq. (|7.ip we obtain 



Combining the above results with Eq. (|2.1ip or 
Eq. (|2.22[) we find the following relations for two strongly 
charged plates with 77 S> 1: 



2T 



WHP 



1 

4' 27 



= L 



2T 



WHP 



2T 



WHP 



a 2 
3 


a 3 
27 


108 J 


= L + 


4 

'/' 


a 2 
3 


a 3 
27 


108 y 


= L + 


2£i(a) 



iplates q 



-plates 2q 



■, —plates 2q : 
(7.8) 

These result hold in the whole strongly charged regime 
where the Gouy-Chapman length Iqc is much shorter 
than the Debye length £db- This includes many exam- 
ples in biological physics. For a given surface charge den- 
sity 77, the a — L curve corresponding to Eqs. (|7.8p can 
be obtained by a rigid shift of the curves for Eqs. (|5.6[ 
5 .91 15. 9p corresponding to infinitely charged plates along 
the L axis. 

Finally for the general case of arbitrary surface charge 
density, we can use Eq. ([7TTj) and Eqs. (|7.2al I7.2b[) to 
numerically solve for a with given L,rj. For example, for 
a given distance between plates L = 2£db, we plot the 
reduced net interaction 5a as a function of the absolute 
value of the surface charge density jTyj in Fig. [S] for all 
three cases. One can explicitly see how the interactions 
saturate as the surface charge density |?7| approaches 
infinity. One can also see how divalent counter-ions 
dramatically reduces the repulsion between two charged 
plates, comparing with the case of monovalent counter- 
ions. With slightly more efforts, we can also plot a as a 
function of separation L for given surface charge density 
77. We shall however not discuss this in detail here. 



VIII. CONCLUSION AND 
ACKNOWLEDGEMENT 

We have discussed in great detail the exact solution of 
the Poisson-Boltzmann equation in the two-plate geom- 
etry, both in q : —q symmetric and in 2q : — q asym- 
metric electrolytes. The Weierstrass elliptic represen- 
tation of the potential has numerous advantages: 1) it 
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Appendix A: Derivation of Poisson-Boltzmann 
Equation 

There are many different derivations of the Poisson- 
Boltzmann equation in the literature. To make this work 
self-contained, we present two different derivations based 
on the variational principle. The equations derived from 
these two approaches are however slightly different. Our 
derivations have the merit of explicitly treating the spa- 
tial variation of dielectric constant and the electrostatic 
boundary conditions. 



Derivation From the Grand Canonical Ensemble 



FIG. 6: The net pressure (a — oto) as a function of the surface 
charge density \r/\ for given plate separation L = 2, obtained 
by solving Eqs. (|7.1l I7.2|l All quantities are dimensionless. 
Red: charged plates in a symmetric electrolyte. Black: pos- 
itively charged plates in a 2q : —q asymmetric electrolyte. 
Green: negatively charged plates in a 2q : —q electrolyte. 
Note how multi- valence counter-ions can dramatically reduce 
the pressure between two plates. 



yields simple analytic results for the one-plate problem 
in 2q : — q asymmetric electrolytes; 2) it yields novel 
and exact asymptotic results in various regimes, both in 
symmetric electrolytes and in 2q : —q asymmetric elec- 
trolytes; 3) it allows us to see explicitly how multi- valence 
counter-ions significantly reduce the repulsion between 
two likely charged plates, as well as how the renormal- 
ized surface charge density depends on the bare surface 
charge density for the one plate problem. 

Our exact solutions of PB may also tell us what to be 
expected beyond the Poisson-Boltzmann theory. There 
are at least three important issues that are ignored by 
the PB. Firstly, when the surface charge density is suffi- 
ciently high, the counter-ion density near the plates may 
become so high that they may crystallize on the plate 
surface. This possibility has been extensively explored 
by many authors, but has not been completely clarified. 
Secondly, if the ions in the electrolyte are of multiple va- 
lence, the Bjcrrum length can be much longer than the 
ion diameter, and most ions form neutral bound pairs. 
The response of these pairs to the external charged ob- 
jects is a serious problem, but has not been explored so 
far. Finally, there may be important short range chemi- 
cal interactions between particular types of ions and the 
charged surfaces. These interactions in principle must 
exist in reality, and can be easily taken into account in 
a field theoretic formalism of the problem. Study of the 
latter two issues will be the major mission of our future 
works. 

The author thanks Leo Radzihovsky, Andy Lau, Hon- 
gru Ma, Erik Luijtcn, Michael Brenner, and Anatoly 
Kolomeisky for helpful discussions on the general field 
of electrolytes. 



Let us start with a symmetric electrolyte with some 
charged dielectric solids fixed in the space. The ions move 
inside the solvent but are not allowed to penetrate into 
the solids. The grand canonical partition function for 
the system of ions can be mapped into the generating 
functional of the sine-Gordon field theory: 



C~ l J Dtp, 



-I3H[ V ] 



(Al) 



H[<p] 



dx 



-e(x)(V<p) 2 - 2Tn{x) cosfiqip 



+ i <t> ipadA, 
I an 



where j3 = 1/T, q is the charge of all ions, while C is 
a normalization constant that is irrelevant to our pur- 
pose. f2 is the spatial region occupied by the dielectric 
solids, while 9f2 is the surface of dielectrics. The func- 
tions e(x), n (x) are given by 



e(x) 



ei, x $ Q 



n (x) = 



°' X ln ■ (A2) 



At the saddle point level, n can be understood as the 
density of each specie of ions, so that the total ion den- 
sity is 2n . cr is the surface charge density prescribed 
externally on the surface of dielectrics, assumed to be 
constant through out this work. Because of the imagi- 
nary boundary term, the saddle point value of the order 
parameter ip (and its average as well) is purely imagi- 
nary. The average of itp has the physical significance of 
the average electrostatic potential. A detailed derivation 
will be published elsewhere. 

Minimization of the action over $ = iip, with the 
boundary terms properly taken into account, leads to 
the Sine-Gordon equation 



-ej A$ + 2nq sinh (3q<& = 0, 

-e s A$ = 0, x i n 



x e fi (A3) 
(A4) 



together with the standard electrostatic boundary condi- 
tion on the dielectric interfaces: 



On 



fm.i 



<9$ 

' dn 



-(7 = 0, 



(A5) 
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where the unit normal vector n goes from the solvent into 
the solid dielectrics. Note that the boundary condition 
Eq. (|A5j) is not satisfied by a generic configuration of the 
field tp, but only by the saddle point configuration. Now 
for aquarius electrolytes ej w 80 > e 8 ~ 1, we can ignore 
the term proportional to e s in Eq. (|A5[) . This leads to 
the Neumann boundary condition: 



Q 



<9$ 

dn 



(A6) 



aa,i 



hence the region inside the dielectric solids is decoupled 
from the region occupied by the electrolyte. We shall 
take this approximation in this work. 



2. Derivation from Canonical Ensemble 

We shall present a slightly different derivation of PB 
from the canonical ensemble. First consider the simple 
case of uncharged solid dielectrics inside a solvent with 
no ion. A unit point charge [l6j at y inside the solvent in- 
duces surface charges on the dielectric boundaries. Inside 
the solvent, the electrostatic potential, i.e. the Green's 
function of the corresponding electrostatic problem, can 
be written as 



G(x,y) 



1 



AtteAx 



y\ 



+ x(x,y), 



(A7) 



where the first term is understood as the potential gener- 
ated by the source charge, while the second term \{x, y) 
is the potential by all the induced surface charges, given 
that there is a unit source charge at y. One can eas- 
ily show that the interaction energy between the point 
charge q and image charges it induces is given by 



W{y) = ^q 2 x{y,y). 



(A8) 



It can be proven that the Green's function satisfies the 
following PDE: 



V ■ e(x)VG{x,y) = 5{x-y). 



(A9) 



Now consider the problem of charged dielectric solids 
inside electrolyte. The total energy of all mobile ions 
(with charges and location n) is given by 

,?*,•) + l;^2q?<f>i{ri) + ^2qi0cxt{n), 

i i 

(A10) 

where the sums are over all the mobile ions, and 

Mv) = x(v,v) (AH) 

is the image charge potential acting on the ion at y, while 
</>ext is the potential produced by all fixed charges on the 
surface of dielectrics. We want to calculate the canonical 
partition function [l7j 



Z = 



Jdn 



(A12) 



where the integrals are over the positions of all mobile 
ions, restricted inside the solvent, i.e. not allowed to 
penetrate into the solids. We consider the following vari- 
ational probability distribution function (pdf): 



N- 



P[{r}]=Hp+(x a )l[p-{yp), (A13) 

where we use coordinates x a and yp for positive and neg- 
ative ions. Since all spatial variables arc independent to 
each other in this pdf, we are ignoring all the spatial cor- 
relations between different ions, a hallmark of all mean 
field theories. 

We calculate the variational free energy 



(A14) 



= / i[d n (P[{r}]H + TP[{r}} log P[{r}]) . 

This variational free energy then provides upper bound to 
the exact free energy, according to Feynman's variational 
principle. 

The variational free energy F mr is more conveniently 
expressed in terms of ion number densities 



n±(r) = N±p±(r). 



(A15) 



After some calculation and ignoring terms smaller by or- 
ders of N' 1 , we find 



Fvar = 7 y ij / 2 K(l)-n_(l)]G(l,2)[n + (2)-n_(2)] 
+ Q 0ext(r)[n+(r) -n_(r)] 



<M r ) [n+(r) +n_(r)], 



+ T [n + (r)logn+(Y) + n_(r) log n_(r)] , (A16) 



where symbols 1, 2 stand for r\,r2- All the integrals are 
restricted outside the dielectrics. Note that the image 
charge potential <fii explicitly contributes to the free en- 
ergy. 

Now define the averaged total potential $(r) by all the 
charges, including the fixed surface charges, the induced 
surface charges, as well as the mobile ions. It is given by 



$(1) = gjG(l,2)[ri + (2)-n_(2)] + Ul), 



where 



0(1,2)^(2)4-^(1), 



p(r) = q[n+(r) - n_(r)] 



(A17) 



(A18) 



is the average mobile charge density. Note that in 
Eq. (|A17|) the total potential is defined in the whole 
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space, including the region inside the dielectrics, even 
though the charge density p{x) is nonvanishing only in 
the electrolyte. Varying the free energy Eq. (|A16|) over 
the number densities n± subject to the constraints of 
fixed total numbers of ±q ions, we find the saddle point 
equation as 



n±(r) = n exp 



(A19) 



where n serves as the Lagrange multiplier fixing the total 
number of ions. As r — > oo, both $ and 4>\ approach 
zero, hence n± — > n. Hence n is bulk ion density for each 
specie. Note how this saddle point equation explicitly 
depends on the image potential </>i(r)! Using Eq. (|A19[) . 
the average charge density p(x) due to all mobile ions 
can be expressed as 



p(r) = q[n + (r) - n_(r)] 
1 



(A20) 



-2ngexp 



(3q z Mr) 



sinh[/3<j$(r)] 



Now let the operator —V • e(x)V acting on Eq. (|A17I) . 
and use Eq. (|A9|) . as well as the fact that 



V ■ e(x)V0 ox t = Pext(z), 



(A21) 



bother us. Even though both equations are derived from 
variational principles, the underlying physics is rather 
different. In the canonical ensemble approach, we di- 
rectly vary the spatial pdf of every ion, while in the grand 
canonical ensemble approach, we vary ip, which is an aux- 
iliary field introduced in a lattice model. A more impor- 
tant difference between two approaches is that Eq. (|A16|) 
is real, while the Sine-Gordon action Eq. (|A3[) is complex. 
As a consequence, the canonical ensemble approach also 
provides an upper bound for the free energy, while the 
grand canonical ensemble approaches does not. Other 
than that, however, it is a priori not so clear which of 
these two equations is a better approximation. Eq. (|A3[) 
nevertheless has the advantage of being simpler. 

Quantitatively the difference between these two MFTs 
becomes relevant only when the exponent —/3q 2 (f>i(r)/2 
in Eq. (|A23[) is large. For the simple case of an ion near 
a flat dielectric boundary with e s <C q w 80, we can esti- 
mate the exponent as — q 2 /8iredT, where d is the distance 
from the ion to the plate. For monovalent ion, this factor 
is becomes larger than unity when d is smaller than half of 
the Bjcrrum lengh, 0.35nm. A few times of this distance 
away from the dielectric boundary, we can safely ignore 
the issue of image charge potential altogether. Eq. (|A3[) 
and Eq. (|A23[) then becomes identical. 



we find 



V ■ e(x)V$(x) = p{x) + p oxt (x), 



(A22) 



where p e xt(x) is the externally prescribed surface charge 
density, which is nonzero only on the dielectric interface. 
Let us now consider the regions inside the electrolyte 
and that inside the dielectrics separately. For the re- 
gion inside the electrolyte, e(x) — e;, while p(x) is given 
by Eq. (|A20|) and p e *t{x) — 0, therefore the equation 
Eq. (|A22j) reduces to 



-e/A$ = — 2nqexp 



sinh[/3g$(r)], x ^ fl, 



(A23) 

which is slightly different from Eq. (|A3|) . For the region 
fi inside the solids, e(a;) = e s , p{x) = p cx t(x) = 0, and 
therefore the equation Eq. (|A22[) reduces to the Poisson 
equation Eq. (|A4[) . Finally, the singular surface charge 
density term in Eq. (|A22[) results in the standard dielec- 
tric boundary condition Eq. (|A5[) . 

Substituting Eq. (|A19j) back into Eq. (|A16|) . we find 
that, up to a trivial additive constant, the variational 
free energy reduces to 



[n + (l)-n_(l)]G(l,2)[n+(2)-n_(2)] 



1 J2 



<r 

2 



(A24) 



That there exists two different versions of the sad- 
dle point equation, Eq. (|A3[) and Eq. (|A23|) should not 



Appendix B: Some General Results at Saddle Point 

Level 



In this section, we derive some general results about 
the PB theory in a symmetric electrolyte. We first define 
three useful length scales (l8j : 



Idb 


= K- 1 = 


toe 


= e/q^a 








q 2 



Q 



2/3nq 2 



Debye screening lengtl^pi) 
Gouy-Chapman length,(B2) 
Bjcrrum length. (B3) 



The following relation between Idb and Xbj shall be use- 
ful below: 



&irt DB \ B jn = 1. 



(B4) 



We further define the dimcnsionless versions of the po- 
tential, spatial coordinates as well as the surface charge 
density: 



V> = Pqtp, 



'2/3n q 2 



<y idb 



V2eTn tec 
We shall also define a dimensionless temperature 



T 



1 



A 



47T- 



Bj 



2nl DB i DB 



(B5) 
(B6) 

(B7) 
(B8) 
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as well as a dimensionlcss Hamiltonian 
H = HT/T. 



(B9) 



The dimensionless temperature T is small for diluted 
electrolytes. For example, in the recent experiment 
on highly deionized water by Tata et al fll| . Iub ~ 
500nm, \sj ~ 0.7nm, and T w 1/57; the system is there- 
fore well in the low temperature limit. The problem of 
dilute electrolytes is therefore simple from the field the- 
ory point of view. 

In terms of all these dimensionless variables, the grand 
canonical partition function Eq. (|A1I) can be written as 



Z = C 



Dipe 



(BIO) 



Hlip] 



dx 



-(Vip) 2 — cosV' 



+ i <t r]ip(r) dA. 
'on 



Varying H over * = iip, and taking into account the 
translational symmetry in the plane of plates, we find 
the saddle point equation 



<9?* + smb.* = 0. 



(Bll) 



The boundary condition Eq. (| A6|) reduces to (again in 
the limit e s /ei — > 0) 



<9* 

dn 



(B12) 



boundary 



The saddle point approximation of the grand potential 
is given by: 

g = TlogZw = -#[-i*]T/T, 

(B13) 



H[-iV] = A dz 



--(<9 z *) 2 -cosh* 



(B14) 



where A is the area of the plates, while *o is the elec- 
trostatic potential on each plate, which should be de- 
termined using the boundary condition Eq. (|B11[) . The 
potential is constant on the plate because of the transla- 
tional symmetry. It is no longer constant for the prob- 
lem of non-planar objects. Rather curiously, the action 
Eq. (|B14|) at the saddle point has a negative gradient 
term. 

Even though we will be primarily concerned with 
the mean field theory, for completeness, we show the 
quadratic part of fluctuation Hamiltonian around the 
saddle point. Let ip = — i* + </>. After some simple 
calculation, we find 



S 2 H[^ *] 



dx 



(V0) 5 



1 



cosh*(r) 2 



(B15) 



It is a quadratic theory with a spatially variable mass 
cosh*(r). Interestingly, the mass becomes exponen- 
tially large where the saddle point potential *(r) is large, 
which makes fluctuations quantitatively unimportant. 



We readily obtain the following first integral for PB 
equation Eq. ()B11|) : 



a = cosh*(z) - -(d z V(z)Y 



(B16) 



Let us show that a > 1 corresponds to the case of two 
likely charged plates, while a < 1 corresponds to the case 
of two oppositely charged plates. If two plates carry the 
same charge density, in the middle of two plates, we have 
d z ^(MP) = due to symmetry, where MP stands for 
"mid point" . Hence 



a = coshf(MP) > 1. 



(B17) 



On the other hand, if two plates carry equal but oppo- 
site charges, then if! (MP) = due to symmetry while 
df(MP) ^ 0, hence 



a = cosh 



l -(d z t>(MP)Y < 1. 



(B18) 



When two plates are infinitely separated, both \E' and 
d z if> vanish in the middle, therefore a = 1, regardless of 
the signs of charge on the plates. 

For the case of two likely charged plates, in no where 
between two plates can the potential \& vanish, for this 
would lead to a < 1, which contradicts our earlier ob- 
servation Eq. (|B17|) . Hence the mean field electrostatic 
potential does not change sign. The so-called "charge 
inversion" does not happen at the saddle point level. 

The total ion density at z is given by 

n tot (z) = n (e^ + e^ 9 *) =2n cosh^(z). (B19) 
Infinite away from the plates, i.e., in the bulk, we have 

f(oo) = 0, ntot(oo) = 2 7i, (B20) 
while in the middle of two plates, we have 

ntat(MP) = 2ncosh*(MP). (B21) 

Hence 

n tot (MP) 



a = cosh* (MP) 



ntot(oo) 



(B22) 



is the ratio between the ion density in the middle of two 
plates to that of the bulk. The fact that f(MP) > 
means that there is a surplus of ion density in the middle 
of two plates. 

To find the interaction between two plates, we further 
rescale the coordiate: 

z^Z = z/L, (B23) 
so that the grand potential Eq. (|B14[) becomes 



TQmf 
AT 



dZ 



- — (dzf>f - L cosh* 



2 77*0, 
(B24) 
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where A is the surface area of both plates. The saddle 
point solution ^ generically depends on L. When we 
vary the above expression over L, however, we find that 
all terms involving 6^ , including the surface term, cancel 
at the saddle point. Therefore the variation of the free 
energy as we change the plate separation L is given by 



at 5Gmf 



T J! dz & dz ^ - L cos *) 

J dz (^(d z V) 2 -cosh*) 



L 

= —a5L, 



(B25) 



where we have transformed back to the usual coordinate z 
in the second equality, and have used Eq. (|B16[) . Restor- 
ing the physical units, then, we find that the force per 
unit area between two charged plates is given by fl9^ 



P 



,- 3 i dg 



MF 



-DB 



A dL 



it 
T 



"3 , 
DB 



T a 



4ire 2 DB \ Bj 



(B26) 



In the first equality, the factor £ DB is needed to restore 
the physical unit for P; in the second equality, used was 
Eq. ([BIB]) : the third, Eq. (TBS) and Eq. ([Hi]) . The last 
equality is useful, because it involves only temperature 
and two fundamental length scales for the electrolyte, as 
well as a dimcnsionlcss parameter a. Using Eq. (IB20[) 
and Eq. (|B22|) . we can also write the pressure as 



P = 2nTa = n tot (MP)T 



(B27) 



Therefore, P is essentially the osmotic pressure of the 
ionic gas at the middle of plates. 

The ionic gas outside the plates also exert pressure on 
them, therefore P calculated above is the not net force 
per unit area between two plates. To obtain the net 
interaction, let us introduce another fictitious plate to 
the left of the left plate with a distance that is much 
greater than L. The electrostatic potential is given by 
the same solution, but with the constant a ~ 1. The 
resulting pressure acting on the plate (pointing to the 
right ) is then given by 2nT. Hence the net pressure 
acting on the left plate is 



P nct = 2nT(a - 1) = 2nT5a, 



(B28) 



which vanishes as the separation L between two plates 
goes to infinity, as it should be. The fact that a > 1 
means the interaction between two likely charged plates is 
always repulsive. Like charge attraction docs not happen 
for plates geometry at the level of mean field theory. The 
result has been proven for more general geometry by Ncu 
some time ago [12(- Finally using Eq. (|B4[) . we can also 
express the net interaction as 



Pn 



TSa 



4nl 2 DB X 



Bj 



(B29) 



Appendix C: 2q : —q Asymmetric Electrolyte 

An outstanding advantage of the Weierstrass function 
representation is that it allows an explicit calculation for 
the case of 2 : 1 asymmetric electrolytes. Assume pos- 
itive ions carry charge of 2q, while negative ions carry 
charge —q. To keep overcall charge neutrality, the bulk 
ion densities are n/2 and n for positive and negative ions 
respectively. The Sine-Gordon field theory for asymmet- 
ric the electrolyte is given by the following: 



Z = c- 1 

H[(p] = J dx 
+ i 



-PH[ V ] 



Dip e 



^e(V^) 2 



(CI) 



Tn I + -e~ 2ifiqtp 



d> ipadA, 
Jan 



For the more general case of kq positive ions, e~ 2% ^ qtp j2 
should be replaced by e~ k ' lf3w /k. 

The Debye length for a 2q : — q electrolyte is given by 



V -DB 



in fiq 2 ' 



(C2) 



with n the bulk density of the — q ions. Wc define the 
Bjerrum length Xgj and the Couy-Chapman length the 
same way as in Sec. [B] The analogue of the identity 
Eq. (|B4f becomes 



\2Kt 2 DB \ B3 i 



1. 



(C3) 



We rescale the field tp, the real space coordinates as well 
as the surface charge density in the following way: 



ip = (3qip, 



'if3nq 2 



V 
f 



x 

= w 

-DB 



= 47T 



Xbj 



"DB 



tec 

1 

3n£ DB 



(C4) 
(C5) 

(C6) 
(C7) 



After the rescaling the partition function can be writ- 
ten as 



Z = 
H = 



D^e- H ' T , 
1 



dx 



iijj 



-2iip 



6 



+ i 



(C8) 
rjtp dA. 



<>v. 



The saddle point is again purely imaginary, i\) = —ity, 
with ^ the electrostatic potential at the mean field level. 
The saddle point equation satisfied by VP is given by 



= 0, 



dn 



= V- 



(C9) 
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Linearization of this equation leads to Eq. (JTT3J), as it 
should be. The first integral of this second order ODE 
can again be easily obtained: 



1 



1 



-2* 



(CIO) 



Consider the case of two equally charged plates. At 
the middle between two plates, d z ^(MP) = 0, hence 



= *(MP) 



-2*(MP) _ n tot(MP) 



3 77 



> 1/2, (Cll) 



where 



»tat (MP) 



(C12) 



is the total ion density at the middle of plates. The lower 
bound a = 1/2 corresponds to the case fy(MP) = 0, 
where two plates are infinitely far away. 

Following a similar strategy as in Sec. Sec. [B] we can 
prove the following result for the pressure between two 
plates: 



P = 3nTa 



Ta 



(C13) 



i.e. it is the osmotic pressure of the ion gas in the middle 
of plates. The net interaction between two plates is then 



given by 



Pnet = 3nT {a - 1/2) 



TSa 



M 2 DB \ Bj 



(C14) 



Since a > 1/2, two likely charged plates repel each other. 

In the bulk, * = d z ^> = 0, hence a = 1/2. Therefore 
the osmotic pressure of the ion gas in a 2q : —q elec- 
trolyte, according to Eq. (|C14[) (with a set to 1/2) is 
given by 



Plq-.-q 



T 



8ire 2 DB \ Bj 



(C15) 



By contrast, the osmotic pressure of the ion gas in a q : 
—q symmetric electrolyte, according to Eq. (|B26|) (with 
a set to unity), is given by 



(C16) 



which is twice of Piq-.-q- This result can be easily ob- 
tained by inspecting the expressions for Debye length in 
both cases, Eq. (TO|) and Eq. (j2~T5j) . 
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